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Abstract 

Single W boson production in electron-positron collisions is studied with the L3 
detector at LEP. The data sample collected at a centre-of-mass energy of y/s = 
188.7 GeV corresponds to an integrated luminosity of 176.4 pb _1 . Events with a 
single energetic lepton or two acoplanar hadronic jets are selected. Within phase- 
space cuts, the total cross-section is measured to be 0.53 ± 0.12 ± 0.03 pb, consis- 
tent with the Standard Model expectation. Including our single W boson results 
obtained at lower y/s, the WW7 gauge couplings k 7 and A 7 are determined to be 
k 7 = 0.93 ± 0.16 ± 0.09 and A 7 = -0.31 t&gj ± 0.13. 
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1 Introduction 



Precise measurements of trilinear gauge boson couplings constitute a crucial test of the Standard 
Model of electroweak interactions PQI2|- The studies of single W productiorj^] 

e+e" -> e + v e W- (1) 

by the LEP experiments [SHU] demonstrated that this process provides one of the best exper- 
imental grounds for precision measurements of the electromagnetic gauge couplings of the W 
boson. The cross-section of process (JTJ) depends only on the k 7 and A 7 gauge coupling parame- 
ters jZj which are related to the magnetic dipole moment, /iw = (e/(2mw)) (1 + k 7 + A 7 ), and 
the electric quadrupole moment, Qw = ( — e / m w) ( K 7 — ^t)j °f the W boson. Any deviation 
from the Standard Model predictions k 7 — 1 and A 7 = would indicate that the W boson has 
an internal structure. 

A specific feature of this reaction is a final state positron produced at very low polar angle 
and therefore not detected. Thus the signature of this process is a single energetic lepton, if 
the W~ boson decays into lepton and anti-neutrino, or two hadronic jets and large transverse 
momentum imbalance in case of hadronic W~ decay. 

In this article we report on the measurement of the cross-section of single W boson pro- 
duction at a centre-of-mass energy of y/s = 188.7 GeV, denoted hereafter as \fs = 189 GeV. 
Combining these results with those on single W boson production obtained at lower -y/i |4 , we 
derive significantly more precise values for the gauge couplings k 7 and A 7 . 



2 Data and Monte Carlo Samples 

The data were collected by the L3 detector [E] at LEP in 1998 with an integrated luminosity 
of 176.4 pb- 1 . 

For signal studies samples of e + e~ — > e + v e ii' events are generated using both the GRC4F [U] 
and the EXCALIBUR JUj Monte Carlo generators. For the background studies the following 
Monte Carlo programs are used: KORALW P] (e+e" -> W+W" -> ffff), KORALZ [T2] 
(e + e~ -> /i + /i"(7), T + T-(i)), BHAGENE3 [13! and BHWIDE [H| for large angle Bhabha 
scattering (e + e~ — > e + e~(7)), TEEGG |T3] for small angle Bhabha scattering (e + e~ — > e + e~7), 
PYTHIA PS] (e+e" -> qq(7)), DIAG36 [TTj and PHOJET [TS] for leptonic and hadronic 
two-photon processes, respectively, and EXCALIBUR and GRC4F for other 4-fermion final 
states. 

The Monte Carlo events are simulated in the L3 detector using the GEANT 3.15 pro- 
gram PH] , which takes into account the effects of energy loss, multiple scattering and showering 
in the detector. The GHEISHA program [201 ^ s use d to simulate hadronic interactions in the 
detector. 



3 Signal Definition 

The signal definition used here is unchanged with respect to our previous publications [HJ[3]. 
The signal is defined as e + e~ — > e + v e ff' events that satisfy the following phase-space cuts: 

^The charge conjugate reactions are understood to be included throughout the paper. 
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I cos 6> c + 1 > 0.997 
min(Et,Ep) > 15 GeV (2) 
| cos6* e - | < 0.75 for e + u e e~u e events only, 

where 9 e + is the polar angle of the outgoing positron, and Ef and Ef, are the fermion ener- 
gies. The final states e + e~ — > e + ^ e ff' that do not satisfy these conditions are considered as 
background; they consist mostly of the reaction e + e~ — > W + W~. In the case of the e + u e e~u e 
final state the additional angular cut reduces contributions from processes where the v e v e pair 
originates from the decay of a Z boson. 

Inside the phase-space region (J2J), single W production dominates since it peaks strongly at 
| cos #c+| ~ 1. On average it accounts for 90% of all events in this region, the remaining 10% 
being mostly non-resonant final states. The purity depends slightly on the flavour of the f r 
pair from W~ decays. For the e + i' e e~u e final state, the purity is 75%. 

For comparison with theory, the cross-sections for this signal definition are calculated with 
a statistical precision from 0.2% to 1.0% using the Monte Carlo generators GRC4F and EX- 
CALIBUR. The main difference between the two generators is the approximation of massless 
fermions used in EXCALIBUR. The reduction of theoretical uncertainties on predictions for 
single W production is subject to ongoing theoretical efforts j2lj. With respect to this dis- 
cussion, we estimate the theoretical uncertainty on the calculated cross-section for single W 
production to be of the order of 7%. This includes the effect of using a smaller electromagnetic 
coupling accounting for the low momentum transfer of the photon in single W production and 
taking into account QED radiative corrections expected for such a t-channel process. Rescaling 
with respect to a smaller electromagnetic coupling constant would lower the cross section by 7% 
to 10%, whereas replacing the energy scale used in the structure function approach to calculate 
QED corrections with the correct physical scale is expected to increase the cross-section by 
about 5%. 

4 Analysis 

All decay modes of the W boson are considered, leading to the following experimental signatures: 
events with two hadronic jets and large transverse momentum imbalance and events with single 
energetic electrons, muons and taus. In the following, efficiencies are quoted with respect to 
the phase-space cuts (j2j), with errors due to Monte Carlo statistics. 

4.1 Hadronic Final States 

The selection of candidates for the hadronic decay W~ — > qq' of single W boson production 
is based on the following requirements: two acoplanar hadronic jets, no leptons, and large 
transverse momentum imbalance. 

High multiplicity hadronic events with at least five charged tracks are selected with an 
energy deposition in the electromagnetic calorimeter greater than 10 GeV and large total visible 
energy, E vis > 60 GeV. All calorimetric clusters in the event are forced to two hadronic jets 
using the DURHAM jet finding algorithm [22j. The invariant mass of the jet-jet system must be 
in the range 40 GeV < M v - ls < 110 GeV. The energy Efb in the forward-backward luminosity 
calorimeters, covering the angular range 0.025 rad < 9 < 0.151 rad, where 9 is measured with 
respect to the incoming electron or positron, is required to be smaller than 65 GeV. These cuts 
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reduce contributions from the purely leptonic final states e + e~ — > e + e~(7), /i + /i~(7), t + t~(j) 
and hadronic two-photon interactions while keeping a significant fraction of hadronic events 
from the processes e + e~ — ► qq(7), e + e~ — > W + W~ and e + e~ — > ZZ. 

To reject events from the two-fermion production process e + e~ — > qq(7), the missing trans- 
verse momentum must exceed 15 GeV. The missing momentum vector must be at least 0.30 rad 
away from the beam axis and the energy in the ±0.22 rad azimuthal sector around its direc- 
tion must be below 10 GeV. In addition, the opening angle between the two jets in the plane 
perpendicular to the beam axis must be smaller than 3.0 rad. 

Events containing identified electrons, muons or photons with energy greater than 20 GeV 
are rejected in order to suppress the remaining background from e + e~ — > W + W~ events where 
one of the W bosons decays into leptons. In order to remove part of the remaining t + u t qq' 
final states with the r lepton decaying hadronically, the event is forced to a three-jet topology 
using the DURHAM algorithm. The solid angle, Q, defined by the directions of these jets is 
required to be smaller than 5.5 sr. 

At 189 GeV centre-of-mass energy, the production of two on-mass-shell Z bosons is a source 
of an additional background. A decay of one Z into neutrinos accompanied by a hadronic decay 
of the other Z leads to an event signature close to that of the signal process ((H). To suppress 
this background, we make use of the fact that the Z bosons produced in pairs are close to the 
kinematic threshold and therefore have low momenta. Thus we ask the velocity of the detected 
hadronic system, (3, calculated as the ratio of the missing momentum to the visible energy to 
be greater than 0.35. 

A total of 216 events are observed in the data, with 35.9 events expected from the signal 
according to the EXCALIBUR Monte Carlo prediction and 179 from the background sources 
distributed as follows: W + W~ production and non resonant e + z/ e qq' final states (166.2 events), 
ZZ production (7 events), two-fermion final state processes (4.2 events) and two-photon in- 
teractions (1.5 events). The contributions to the background from other processes are found 
to be negligible. The number of observed events is in good agreement with the expectation. 
The signal selection efficiency is determined to be (50.5 ±0.7)% using the EXCALIBUR Monte 
Carlo and (51.8 ± 1.3)% using GRC4F Monte Carlo. 

In order to differentiate further between the signal and the W + W~ background, a discrim- 
inant variable NN out is constructed using a neural network. Nine variables are combined in 
a feed-forward neural network [23], with one hidden layer and one output node. The input 
to the neural net includes the spherocity of the event, the invariant mass of the two jets, the 
masses of the two jets, the velocity j3 of the hadronic system, the solid angle Q, the resolution 
parameters 2/23 and 2/34 of the JADE jet finding algorithm [21] for which the number of jets in 
the event changes from two to three and three to four, respectively, and the ratio of the mass 
to the energy of the least energetic jet after forcing the event to three jets. As an example, the 
distribution of the solid angle Q for the selected events is shown in Figure Q^a). The use of the 
neural network increases the signal fraction in the selected sample to approximately 60% for 
large neural network output values, as shown in Figure ^b). 

The cross-section of the process (1) for hadronic decays of the W boson is determined by 
a binned log-likelihood fit to the neural network output distribution of Figure lU^b), assuming 
Poisson statistics. The background contributions are fixed to the corresponding Standard Model 
Monte Carlo predictions. The fitted cross-section of the hadron channel is found to be: 

cr(e + e" -> eu e qq) = 0.41 ± 0.11 pb . 

A similar result is obtained if the GRC4F Monte Carlo is used for the simulation of the signal. 
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This result has to be compared with the expected signal cross-section of 0.40 pb predicted by 
EXCALIBUR and 0.38 pb predicted by GRC4F. 

As a check of the analysis procedure, a fit of the total W + W~ cross-section, keeping the single 
W contribution fixed to the EXCALIBUR Monte Carlo prediction, gives the value 16.9 ± 1.5 pb 
in good agreement with the Standard Model expectation of 16.7 pb predicted by KORALW. The 
result of the fit with both W + W~ and single W production cross-sections as free parameters is 
in good agreement with the results above. 

4.2 Leptonic Final States 

The distinct feature of the process e + e~ — > e + u e W~ , W~ — > t~i>i is a high energy lepton from 
W decay with no other significant activity in the detector. 

Events with one charged lepton (electron, muon or tau) with an energy of at least 15 GeV 
are selected. The lepton identification is based on the energy distribution in the electromagnetic 
and hadron calorimeters associated with the trajectory of charged particles. The total energy, 
^vis, is calculated as the sum of the lepton energy, E e , measured as discussed below, and the 
energies of all neutral calorimetric clusters in the event. No other charged particle activity is 
allowed. The ratio E^/E v - m is required to exceed 0.92 in order to suppress background from 
two-fermion production e + e~ — > £ + £~(7). In addition, the energy in the ±0.22 rad azimuthal 
angle sector along the missing energy direction must be below 1 GeV (0.2 GeV for muons). The 
energy in the forward-backward luminosity calorimeters, Ep-Q, must be less than 70 GeV. The 
accepted background is dominated by two-fermion production processes, especially radiative 
Bhabha scattering in the case of the single electron final state, e + v e e~v e . Moreover, significant 
contributions are due to 4-fermion final states that include two neutrinos and fall outside the 
signal definition (J2J). 

4.2.1 Single Electron Final States 

Electrons are identified as clusters in the electromagnetic calorimeter consistent with an elec- 
tromagnetic shower shape matched in azimuthal angle with a track reconstructed in the central 
tracker. For the single electron final states, the electron energy, E e , measured in the electro- 
magnetic calorimeter, must exceed 20 GeV and the polar angle 8 e must satisfy the condition 
| cos 8 C \ < 0.7. This requirement reduces the contribution from Bhabha and Compton scattering 
and from the process e + e~ — > e + e~i > i> where the e + e~ pair originates from a low- mass virtual 
photon. The e + e~ — > e + e~(7) events constitute the dominant contribution to the selected sam- 
ple. The requirements E e /E V i S > 0.92 and Ep-g < 70 GeV reduce significantly this contribution. 
The acoplanarity angle between the direction of the electron and the most energetic neutral 
cluster (if any) must be greater than 0.14 rad. 

A total of 8 events are observed with 11.9 expected from the Standard Model including 7.9 
events expected from the signal as predicted by the EXCALIBUR Monte Carlo. The energy 
spectrum of the selected events is presented in Figure a). 

The trigger efficiency is found to be 92% from a control data sample. The signal selection 
efficiency is estimated to be (80 ±2)% using the EXCALIBUR Monte Carlo program and 
(78 ±3)% using GRC4F. The major sources of the efficiency loss are due to the requirements 
| cos#o| < 0.7 and E c > 20 GeV. A binned log-likelihood fit to the electron energy spectrum 
results in: 

a(e + e" -»■ eu e eu e ) = 0.021±g;g|g pb , 
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to be compared with the signal cross-section of 0.056 pb predicted by EXCALIBUR and 0.054 pb 
predicted by GRC4F. 

4.2.2 Single Muon Final States 

Single muon final states are selected as events containing one isolated muon, identified as a 
reconstructed track in the muon chambers and the central tracker. The muon energy, E^, 
measured in the muon chambers and in the central tracker, is required to exceed 15 GeV. The 
fiducial volume for this analysis is defined to be | cos# M | < 0.85. This latter requirement avoids 
significant decrease of the trigger efficiency for single muons in the angular region | cos # M | > 
0.85. Within the acceptance, the trigger efficiency exceeds 96%. The total calorimetric energy 
not associated to the muon must not exceed 3 GeV. The rejection of background from cosmic 
muons is based on the radial distance of closest approach of the muon to the beam line and a 
match in azimuthal angle of the muon chamber track with a track reconstructed in the central 
tracker. 

A total of 10 events are observed with 9.8 expected from Standard Model processes including 
7.5 events from the signal. The energy spectrum of the selected single muon candidates is 
presented in Figure I2fb). A signal efficiency of (61 ± 2)% is estimated using the EXCALIBUR 
Monte Carlo program and (56 ±2)% using GRC4F. The main source of the efficiency loss is 
due to the geometrical acceptance of the muon chambers. A binned log-likelihood fit to the 
muon energy spectrum results in: 

a(e + e" eu e ^) = 0.0701^27 P b • 

The expected signal cross-section is 0.060 pb according to EXCALIBUR and 0.059 pb according 
to GRC4F. 

4.2.3 Single Tau Final States 

Single tau final states are selected as events containing one low-multiplicity hadronic jet. The 
calorimetric energy associated with the r jet, E T , must exceed 15 GeV. The number of tracks 
reconstructed in the central tracking system must be either 1 or 3. 

A total of 4 events are observed with 3.6 expected from the Standard Model processes 
including 2.1 events from the signal. The energy spectrum of selected single tau candidates is 
presented in Figure|2[c). The signal efficiency is estimated to be (26 ±2)% using EXCALIBUR 
and (29 ±2)% using GRC4F. The trigger efficiency was studied and found to be in excess of 
98%. A binned log-likelihood fit to the tau energy spectrum yields: 

a(e + e" ez/ e rz/ r ) = 0.066^0 pb . 

The predictions for the signal cross-section are 0.060 pb and 0.059 pb according to EXCALIBUR 
and GRC4F, respectively. 

5 Systematic Uncertainties 

In case of the hadronic decay of the single W, the differences of the EXCALIBUR and GRC4F 
signal modelling are taken into account in the systematic uncertainty of the cross-section mea- 
surement. This systematic uncertainty is found to be approximately 3%. In addition, the 
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parameters describing the neural network structure are varied and the analysis is repeated to 
allow an estimation of the uncertainty due to the choice of the network, yielding a contribution 
of 3%. Compared to this, detector effects, studied by smearing and shifting the kinematic 
variables that are fed into the network within the experimental resolution, have a negligible 
impact on the result. 

In the lepton channel, the dominant systematic uncertainty of approximately 4% arises 
from the signal modelling comparing the signal efficiencies estimated using EXCALIBUR and 
GRC4F. The uncertainty due to the identification of leptons is studied using control data 
samples of two-fermion production and is found to be less than 1.5%. 

The uncertainty due to the Monte Carlo signal statistics ranges from 1% to 3% on the 
cross-section depending on the decay channel. The systematic uncertainty on the expected 
number of background events is essentially due to the limited Monte Carlo statistics, with 
smaller contributions from the uncertainties on the cross-sections and the selection efficiencies 
for the background processes. The overall uncertainty on the total number of background 
events ranges from 3% to 4% in the individual channels. These uncertainties are uncorrelated 
among individual channels and different centre-of-mass energies and have negligible impact on 
the final results. Taking into account all the contributions, the systematic uncertainty of the 
cross-section measurement amounts to 5% for the hadronic decay channel and 6% overall. 

6 Results 

6.1 Total Cross-Section 

The total cross-section of single W production is determined from a binned likelihood fit to 
the distributions of the neural network output presented in Figure Hfb) and the lepton energy 
spectra shown in Figures a) - (c). The sum of the different lepton energy distributions is 
presented in Figure I2fd). The background shapes and normalisations are fixed to the Monte 
Carlo prediction. The fitted signal cross-section, cr(e + e~ — > ez/ e W), corresponds to that of 
the process e + e~ — > eu e ii, where ff denotes a sum of all ivi* and qq final states satisfying the 
phase-space conditions (2). The total single W boson cross-section at \/s = 189 GeV is then 
determined to be: 

a(e + e" -> ez/ e W) = 0.53 ±0.12 ± 0.03 pb , 

where the first error is statistical and the second systematic. The measured cross-section value 
is consistent with the Standard Model prediction of 0.57 pb calculated with EXCALIBUR 
and 0.56 pb calculated with GRC4F. The dependence of the cross-section on the centre-of- 
mass energy agrees well with the Monte Carlo predictions as shown in Figure El including our 
previous measurements at centre-of-mass energies between 130 GeV and 183 GeV [I]. 

6.2 WW7 Gauge Couplings 

The electromagnetic gauge couplings k 7 and A 7 describing the WW7 vertex are determined 
from a binned maximum-likelihood fit similar to the one used for the cross section determina- 
tion. In the fit each Monte Carlo event is assigned a weight that depends on the generated 
4-fermion event kinematics and the values of the gauge couplings k 7 and X T The dependence 
of the background from W + W~ production on the gauge couplings is also taken into account. 
The weight is calculated using the matrix element as implemented in EXCALIBUR, imposing 
constraints on the triple gauge boson couplings k% and Az arising from the SU (2) x U(l) gauge 
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invariance: Kz = gf — tan 2 # w (ft 7 — 1) and Az = A 7 . These constraints affect only the back- 
ground contributions, as the signal process depends on A 7 and k 7 only. The general analysis 
of the remaining C- and P-conserving couplings, gf, A 7 and k 7 , can only be done with full 
consideration of the W + W~ production. In the present analysis we fix the weak charge of the 
W bosons to its Standard Model value, gf — 1, and focus on the electromagnetic properties of 
W bosons. 

The dependence of the coupling determination on the total cross-section for single W boson 
production is tested repeating the likelihood fit with a ±7% variation of the signal cross-section. 
The corresponding systematic uncertainty is ±0.04 for k 7 and ±0.02 for A 7 . Comparing the 
signal description of the two Monte Carlo Generators, EXCALIBUR and GRC4F, shows an 
additional systematic uncertainty of ±0.05 and ±0.06 on k 7 and A 7 , respectively. The agreement 
between the two generators for various anomalous couplings is checked. No coupling-dependence 
of the ratio of the two cross-section predictions was found. 

In addition, the cross-section of the background contributions coming from W + W~ and ZZ 
production are varied to allow an estimation of the systematic uncertainty. The influence of 
these uncertainties on the coupling determination is found to be less than ±0.02. 

The estimated systematic uncertainty is assumed to be Gaussian and fully correlated be- 
tween individual channels. The systematic uncertainty in the cross-section determination is 
taken into account by convolution of the likelihood function with a Gaussian in the fit. 

For the fit to the couplings, we combine the single W data at yfs = 189 GeV presented 
here with our single W data already published [I]. Single W boson production is particularly 
sensitive to the gauge coupling k 7 . Thus, this coupling is determined in a fit fixing A 7 = to 
be: 

k 7 = +0.96l{j;^ ±0.09, 

where the first error is statistical and the second systematic. Fixing k 7 — 1 and performing a 
fit for A 7 yields: 

A 7 = -0.26lg;^ ±0.13. 
Varying both couplings A 7 and k 1 freely in the fit yields: 

k 7 = ±0.93±0.16±0.09 
A 7 = — 0.31±8:f| ±0.13 , 

with a correlation coefficient of ±37%. The 68% and 95% confidence level contours on k 7 and A 7 
are shown in Figure 0] The results are consistent with the absence of anomalous contributions 
to WW7 couplings. The limits on k 7 and A 7 at 95% confidence level are: 

0.56 < k 7 < 1.29 for A 7 = 
-0.67 < A 7 < 0.59 for k 7 = 1 . 

These results represent a major improvement in the accuracy on the triple gauge boson couplings 
k 7 and A 7 compared to our previous publications on single W boson production 0111 • They 
are complementary to measurements based on W + W~ production at LEP |6ll2"5] or determined 
at the Tevatron pp collider [2E] . 
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Figure 1: Distribution of (a) the 3-jet solid angle fl and (b) the neural network output for the 
selected hadronic events. The points are data, the hatched histograms represent the background 
and the open histograms show the expected signal from ez/ e qq final states as predicted by the 
EXCALIBUR Monte Carlo. 14 



L3 

(a) 




• Data 1 89 GeV 

□ ev e ev e 

C3 Background 



> 

CD 

CD 
o 



co 

H — ' 

CD 
> 
LU 



20 40 60 80 100 

Electron Energy [GeV] 



L3 

(b) 



• Data 189 GeV 

□ ev e ^ 

E3 Background 



o o 



20 40 60 80 100 

Muon Energy [GeV] 



L3 

■(c) 



• Data 1 89 GeV 

□ ev e xv T 

C3 Background 



20 40 60 80 

Tau Energy [GeV] 



100 



> 

CD 

CD 
o 



CO 
-t— ' 

CD 
> 
LU 



12 



10 



2 




L3 

(d) 




• Data 189 GeV 

□ ev e lv, 

E3 Background 



20 40 60 80 100 

Lepton Energy [GeV] 



Figure 2: The energy spectra of the selected (a) single electron, (b) single muon, and (c) 
single tau candidates. The sum of the different lepton energy spectra is shown in (d). The 
points are data, the hatched histograms correspond to the background contribution. The open 
histograms show the expected signal from ev e £i>e final states as predicted by the EXCALIBUR 
Monte Carlo. 
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Figure 3: The measured cross-section of single W production within our phase-space cuts 
as a function of the centre-of-mass energy. The solid and dashed lines show predictions of 
the GRC4F and EXCALIBUR Monte Carlo programs, respectively. The estimated theoretical 
uncertainty of ±7% is indicated by the band. 
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Figure 4: The contours corresponding to 68% and 95% confidence level in the « 7 — A 7 plane. 
The point indicates the global minimum from the 2-parameter fit, to be compared with the 
Standard Model prediction indicated by the star. 
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